Q 
nJ 

hi 



U 



GET P. N 



= 0 



\0 



1 = 1 + 1 

G = GCD ( P, THOUSAND_PRIMES [I]); 



<^G>1 



110 



RETURN FALSE 



1 




RETURN TRUE 





GETU.V 



GET P. M, E 



DETERMINE M'^ MOD P 
BY CALCULATING 
GCD (P. M) 



<5 



APPLY SIGNED DIGIT 
ALGORITHM TO E 



|6 



CALCULATE M^ MODP 



Fig i> 



U3 = U;V3 = V; 
U2 = 0:V2 = 1: 





RETURN 




U3, U2 



Q = LU3A/3j; 
T3 = U3-Q*V3; 
T2 = U2-Q*V2; 
U3 = V3;V3 = T3; 
U2 = V2;V2 = T2; 




GETU.V&M 



U3 



= U;V3 = V; W'*^ 



D = 8 = 32 ]>^ | D = 5 



[A,B] = TOP2M(U3,V3,M); 




▼ y 



LOG(U3)-LOG(V3)>D 




Q1 = LA/BJ; 
T3 = U3-Q1 *V3 



[A, B] = TOP2M(|T3|. V3. M); <-<S7v3> 
Q2 = Q1- 1a/BJ+1); — ^ 



calculate 

XA + YB<2^^ 



r 



Q2 = Q1 - 1 



T3 = U3-Q2*V3 



N 



T3 = XU3 + YV3; 
U3 = V3; V3 = |T3| 



Q = LU3/V3j 
T3 = U3-Q*V3 



U3 = V3;V3 = T3; 





^^^^ |y 



Y = GCD (U3,UM0D U3); 
Y = GCD (Y.VMODY); 



RETURN Y 



GETA,B,M ^ 


1 


r 




M = 32 




V 



[X, Y] = MOD_EUCLID_GCD(A, B) 



[X. Y] = MOD_LEHMER_GCD{A, B) 



RETURN X.Y 



riG. 5" 



GETU3.V3,M 



4^^ 



PB1 = LOG (U3) 



BP2 = BP1 -2M + 1 



RETURN U3[BP1:BP2], V3[BP1:BP2] 



FIG, 4A 



GETU.V -6?/0 



U3 = U; V3 = V; 
U2 = 0;V2 = 1; 




L0G(U3)-L0G (V3)>32 




[SU, SV] = TOP2M(U3, V3, 16); 
D = LOG (U3)-32; 
J = 0; DONE = FALSE; 
X[0]=1;X[1] = 0;Y[0] = O.Y[1] = 1; 



Uf2A 



Q = LSU/SVj; 
X[J+2] = X[J]-X[J+1]; 
Y[J+2] = Y[J]-Y[J+1]; 
T = SU.Q*SV; SU = SV; SV = T; 
J = J+1; 



D0NE = SV<-Y[J+1]0R 
SU-SV<X[J+1]-X[J]0R 
SU < 



J EVEN > ^ 



r 



D0NE = SU<-X[J+1]0R 
SU-SV<Y[J+1]-Y[J]0R 
SU<2D; 



< DONE = TRUEf > 



T3 = X(J-1] * U3 + Y[J-1] * V3; T2 = X[J-1] * U2 + YO-1] * V2; 
V3 = X[J] * U3 + Y[J] * V3; V2 = X[J] * U2 + Y[J] * V2; 
U3 = T3; U2 = T2; 



Q = LU3/V3j; 
T3 = U3-Q*V3; T2 = U2-Q*V2 
U3 = V3; V3 = T3; U2 = V2; V2 = T2 



I 



PIG. [s> 



<S>ry 



U1=(U3-U2*V)/U; 
RETURN U1.U2 



GETU.V 



7 / O 



\ 



U3 = U;V3 = V: 
U2 = 0;V2 = 1; 



•i 



U3 < 



7 li? 




i 



U1=(U3-U2*V)/U: 
RETURN U1.U2 



Q = LU3/V3j; 
T3 = U3-Q*V3 
T2 = U2-Q*V2; 
U3 = V3;V3 = T3i 
U2 = V2;V2 = T2; 



7/6 



m 



GETU,V&M 



U3 = U;V3 = V; ^^/^^ 
U2 = 0; V2 = 1; 



[A, B] = TOP2M(U3, V3, M); 




LOG(U3)-LOG(V3)>D 




Y 





Q1 = LA/BJ; 


1 — ► 


T3 = U3-Q1*V3 







r 







CALCULATE 
XA + YB<2^^ 



^5? 



1. 



Q = Q1 



■( [t3|>V3^ 



[A, B] = T0P2M(|T3|,V3. M); 
Q2 = Q1- (LA/BJ+1) 



Q2 = Q1-1 



T3 = U3-Q2*V3V ^^ 



T3 = XU3 + YV3; 
T2 = XU2 + YV2: 



T2 = -T2 <- <T3 < 0> 



Q-Q2 U^:^<ii7oS>-^ Q = LU3/V3j 
Q-Q2 ^--^1310;^^ T3 = u3.Q.v3 



U3 = V3; V3 = |T3| 
U2 = V2;V2 = T2 



12 = 


U2-Q*V2; 


U3 = 


V3;V3 = T3; 


U2 = 


V2;V2 = T2; 




V3 = 0 



Y = GCD (U3, U MOD U3); X = LU3A'J; 
COMPUTE X-1 WRT U BY GCD (U. X); 
U2 = (U2*X-i) MODU 
Y = GCD(Y.VMODY); 



